HD-THEP-04-44 



Cosmological bounds on the equation of state of dark matter 

Christian M. Miiller 1 

1 Institut fur Theoretische Physik, Philosophenweg 16, 69120 Heidelberg, Germany 

(Dated: October 26th 2004) 

In this exploratory study, we investigate the bounds on the equation of state of dark matter. 
Modeling dark matter as a fluid component, we take into account both positive and negative fixed 
equations of state. Using CMB, supernovae la and large scale structure data we find constraints on 
the equation of state in a modified ACDM cosmology. We obtain —1.50 x 10 -6 < iuj m < 1.13 x 10 -6 
if the dark matter produces no entropy and —8.78 x 10 -3 < Wdm < 1-86 x 10 -3 if the adiabatic 
sound speed vanishes, both at 3 a confidence level. 
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I. INTRODUCTION 

In this note, we would like to present bounds on the 
dark matter equation of state derived from cosmological 
observations. In other words: given the information from 
present cosmological observations, how cold is cold dark 
matter? In this context, we do not want to limit our- 
selves to a non-negative equation of state since the na- 
ture of dark matter is not yet clear. There are a number 
of particle dark matter candidates and numerous exper- 
iments attempting to detect dark matter particles ( see 
and references therein). However, dark matter may 
not be a particle at all. We will try to keep an open 
mind and check if a negative equation of state can be 
ruled out by cosmological observations. Previous stud- 
ies focused on the power spectrum properties of warm 
dark matter [j| and mixed dark matter models Q, but 
we would like to emphasize that the present work does 
not consider the bounds on warm dark matter. We will 
consider both positive and negative equations of state 
Wdm and show that the equation of state of dark matter 
is already strongly constrained by current observations 
of the CMB, supernovae la and large scale structure. We 
will not concern ourselves with bounds from other than 
cosmological observations. 

This is meant as an exploratory study, not as rigid 
modeling of dark matter, and we will therefore work only 
from a fluid perspective and leave the question open how 
one could obtain a negative equation of state of dark mat- 
ter. For clarity, we will only allow for a constant equation 
of state Wdm , even though from a particle perspective it is 
clear that Wdm varies with time for warm and mixed dark 
matter models. We have chosen this approach as there 
is no simple possibility to give a corresponding particle 
motivation for Wd m < 0, as specifying the time evolution 
of Wdm would be model dependent. 

In models where dark matter interacts with other com- 
ponents of the universe, such as in coupled quintessence 
models one may obtain a negative equation of state 
for dark matter. It may be possible to obtain Wd m < 
by other methods as well, but we are not aware of any 
such model. 

In this work we will use two simple models for dark 



matter: one with no entropy production and one with 
vanishing adiabatic sound speed, both with fixed equa- 
tion of state. 

We obtain bounds for a constant equation of state of 
dark matter of-1.50xl0~ 6 < Wdm < 1.13 x 10 -6 if there 
is no entropy production and —8.78 x 10~ 3 < Wd m < 
1.86 x 10~ 3 if the adiabatic sound speed vanishes, both 
at 3 a confidence level. 

For this investigation we will assume that the universe 
is flat and contains a cosmological constant type dark 
energy component with equation of state Wde = — 1, dark 
matter with a variable equation of state Wdm, baryons, 
photons and massless neutrinos. We do not include the 
tensor part in our analysis. In the conclusions, we will 
address the issue what will change if we relax the flatness 
assumption and if we have an equation of state different 
from —1. 

The plan of this note is as follows: in Section [HI we 
will introduce the perturbation equations for dark mat- 
ter with an arbitrary equation of state (which we will 
from now on call modified dark matter for brevity). In 
Section ITTT1 we introduce a model with vanishing entropy 
production, in Section llVI we give a different model with 
vanishing adiabatic sound speed. In Section[V]we present 
the bounds on Wdm obtained from a Markov Chain Monte 
Carlo simulation, followed by our conclusions. 



II. PERTURBATION EQUATIONS 

In the following, we will use the notation of 5]. Since 
this reference may not be readily available, we will shortly 
introduce our notation. The full energy momentum ten- 
sor for a fluid with equation of state w may expressed 
by (we are suppressing the species index for notational 
convenience) 
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where the Q(k,x) are eigenfunctions of the Laplace- 
Operator, V 2 Qk{x) = —k 2 Qk{x) and in spatially flat 
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universes Q = exp(ikx). The (gauge dependent) vari- 
ables S, v, B, 7T£, and IT are defined by these equations. 
These quantities may then be combined to form the gauge 
independent quantities A g , the density perturbation on 
hypersurfaces of constant curvature perturbation, V, the 
velocity and II, the anisotropic stress component. We 
will construct a gauge invariant form of the pressure per- 
turbation 7Tl later. 

The perturbation equations for the dark matter com- 
ponent expressed in gauge-invariant variables are then 



5- 



3(c 2 



S a A 
W)~Ag 

a 



+ kV(l +w) + 3-wT = 0, 



V = 



a 

elk 



l)V + fc(# - 3^$) 



wk 



T,a (Pa+Pa)V a 



2M z p k 



-n 



(5) 

(6) 
(7) 



where the sum is over all present species. We will as- 
sume that the anisotropic stress vanishes for dark mat- 
ter, Udm = 0. $ and 'i? are the gravitational potentials 
where $ = — ^> — J2 a il a Il a . The sound speed is given 
by c 2 = p/p. Note that this is not the adiabatic sound 
speed, which is defined by c 2 d — Sp/Sp. T is the entropy 
production rate and is given by 
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This may also be expressed as the difference between 
"background" sound speed and adiabatic sound speed, 
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It will be useful to formulate T in terms of gauge- invariant 
variables. From 5j it may be verified that 



ttl := 7T L + 3— (1 + u>)-£;"V 
w a 
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is the gauge-invariant pressure perturbation. Hence 

r = tt z -^[A_- 3(1 +«;)*]. (11) 



III. DARK MATTER WITH NO ENTROPY 
PRODUCTION 

In this section, we will assume that T = 0. For a con- 
stant equation of state of dark matter, the sound speed 



1 Please note that w is the equation of state of dark matter and 
that the species indices are suppressed. 



40 



31) 



< 20 



10 



— w=0 
— w=-0.01 

w=0.04 (x 10) 



FIG. 1: Evolution of the energy density perturbation of dark 
matter A 9 in the sub-horizon regime for several equations of 
state in a universe containing only dark matter with V = . In 
the case of w = 0, the energy density perturbation grows oc r 2 
(black, straight). For < w < 1/3 (blue, dotted) one obtains 
decaying oscillations while for w < the density perturbation 
grows rapidly (red, dashed). 



is given by 



cl = P/p = w. 



(12) 



If we allow for negative equation of state the square of 
the background sound speed may thus become negative. 
There is of course a question whether or not this assump- 
tion is reasonable, but given the unknown nature of dark 
matter we may consider this case and see how well con- 
strained w is. Before we discuss the numerical solutions 
of the perturbation equations, it will be helpful to con- 
sider the solutions for a universe filled only with dark 



matter in the sub- horizon limit, k ^> (- 
©-d ma Y be combined to eleminate V: 
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(13) 



— — w 
2 2 



a . 



3 2 U 



(14) 



The background solution for a universe filled only with a 
modified dark matter component yields 
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For w — the solution of (|13fl and (|14|) in the sub-horizon 
limit is A g = a(r) = t 2 and $ = const., as is well known. 
For the superhorizon regime, A g = const. The solutions 
to these equations in the sub-horizon limit if w ^ are 
plotted in Fig. □ It can be seen that for < w < 1/3 
we obtain decaying oscillations for A g while for w = 1/3 
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FIG. 2: CMB and matter power spectra for a ACDM model 
with modified dark matter, T = 0. The power on small scales 
is much larger than in the w — case if the equation of state is 
negative (red, dashed). One can see the decaying oscillations 
for the case of w > (blue, dotted). As expected, there is 
no difference in the power spectra at very large scales. The 
impact on the CMB spectrum is not noticable at this level. 



one gets oscillations with constant amplitude. For w < 
we find that A g increases rapidly. This is obviously due 
to the negative sound speed, which leads to a growing 
gravitational potential. We can thus already see that 
the equation of state will be strongly constrained in the 
regime w < since we do not observe this excess of power 
on small scales. Also, the growing gravitational potential 
would lead to excessive lensing on small scales, which 
is not observed. At super-horizon scales, A s = const., 
regardless of the equation of state. 

We have computed the resulting CMB andpower spec- 
trum in Fig. |3 with the Cmbeasy software (which is 
based on Cmbfast using the parameters of a best fit 
ACDM model: fl b h 2 = 0.230, fl m h 2 = 0.144, n s = 0.961, 
t = 0.114, h = 0.69. 

As one can verify the impact of varying w on the power 
spectrum is very marked on small scales. We therefore 
expect to obtain tight constraints on w. 



IV. DARK MATTER WITH VANISHING 
ADIABATIC SOUND SPEED 

Given the results of the last section and the problem- 
atic assumption of negative sound speed, we may choose 
a different approach. In fact, T measures the "difference" 
between adiabatic sound speed c 2 d — 8p/Sp and "back- 
ground" sound speed c 2 = pj p. In the previous example, 
we enforced c 2 ad = c 2 by the requirement T = 0. A differ- 
ent requirement would be that the adiabatic sound speed 



vanishes. 

The first problem that arises is that c 2 d is not gauge- 
invariant. We must therefore first specify what we mean 
by vanishing adiabatic sound speed. We may choose a 
hypersurface such that c 2 d = has a definite meaning. 
Since c 2 d = implies Sp = and ttl = Sp/p this leads to 
the requirement that tt^ = on a certain set of hyper- 
surfaces. For simplicity, we choose the Newtonian slicing, 
giving the shear free hypersurfaces, a g = 0. We therefore 
obtain 7r (™ e ™*) — q, which with the definition Eq. 110|) 
gives 

ttl = 0. (16) 

This is true in any gauge because ttl is gauge-invariant. 
Hence choosing 

r = 3(l + w)$-A s , (17) 

we enforce that the adiabatic sound speed vanishes on 
hypersurfaces of isotropic expansion rate. Of course, this 
choice is by no means preferred over any other choice of 
hypersurfaces; we have chosen this one merely for com- 
putational simplicity. 

Solving the perturbation equations for a universe filled 
only with modified dark matter, we obtain evolution of 
super-horizon modes if w ^ 0. There is no exponential 
growth for sub-horizon modes if w < 0. We may conclude 
that this model is well-behaved compared to the T = 
case. 

We have plotted the power and CMB spectra in Fig. [2 
model parameters as for the T = case). As expected, 
the modification has a huge impact on the growth be- 
haviour of fluctuations. 



V. BOUNDS ON THE EQUATION OF STATE 
OF DARK MATTER 

In order to quantify the bounds on the equation of 
state of dark matter, we ran a Markov Chain Monte Carlo 
(MCMC) simulation for T = and c 2 ad = with the An- 
alyzeThis! packa ge M using WMAP TT and TE spec- 
tra [H 113, VSA llJcBI |13| and ACBAR data 
up to I = 2000 as well as the SDSS power spectrum 
(all points with k/h < 0.15 Mpc -1 ) and the SNe la data 
of Riess et al. [l6|- Each run contained ~ 50, 000 points 
after burn-in removal. The model used was a ACDM cos- 
mology with modified dark matter and parameter priors 
as shown in Table [i] 

The resulting one-dimensional marginalized likelihoods 
are displayed in Figure^ The confidence intervals for the 
equation of state are displayed in Table ITU 

The equation of state is quite strongly constrained if 
r = 0, at a level of 10~ 6 . What is somewhat surprising 
is that the likelihood is centered not on w = but on a 
slightly negative equation of state. This may be traced to 
the fact that the SDSS data set we used does not encom- 
pass very small scales and therefore the observations are 



4 




0,01 0,1 
k/h [Mpc '] 



FIG. 3: CMB and power spectrum for a ACDM model with 
modified dark matter, c 2 d = 0. For a positive equation of 
state (blue, dotted) the peak positions of the CMB are shifted 
to smaller scales and have less power than in the w = case 
(black, straight). For a negative w (red, dashed), the peak 
positions are shifted to smaller I. Note also the different peak 
ratios in these cases. The shapes of the matter power spec- 
trum are different in each case, but in contrast to F = 0, there 
is no dramatic difference at small scales with respect to the 
w — spectrum. The power spectra at very large scales are 
different with respect to each other, indicating evolution of 
super-horizon sized modes. 



TABLE I: Flat priors used for the MCMC simulations 



parameter 


min 


max 


Q. b h 2 


0.016 


0.03 


n m h 2 


0.05 


0.3 


h 


0.60 


0.85 


n s 


0.8 


1.2 


T 





0.9 


w dm for T = 


-4 x 10~ 6 


5 x 10" 


w dm for <? ad = 


-0.02 


0.02 



(i 



blind to the very strong increase in power at very small 
scales if w < 0. The constraints for the c^ d = case are 
less stringent, in line with our expectations. Here, w is 
only constrained at a level of 10 -3 . 



VI. CONCLUSION 

Since this is an exploratory study we have not at- 
tempted to formulate a model with realistic variation in 
the equation of state but chose to show the main effects 



TABLE II: Confidence intervals 





68.3 % 


95.4 % 


99.7% 


w dm for r = (xHT 6 ) 


—0.114 
-0.929 


+0.462 
-1.19 


+1.133 

-1.50 


w dm for c 2 ad = (xlO" 3 ) 


+0.098 
-3.59 


+0.890 
-6.58 


+ 1.86 
-8.78 


and bounds derived from 


current 


observations 


using a 



constant equation of state. For the r = model, it is 
clear that the main constraint comes from the matter 
power spectrum. If we would include measurements at 
smaller scales, the constraints on the equation of state 
would be even more restrictive. We may conclude that 
the r = model is unlikely to be realistic. The situa- 
tion is not so clear for the c ad = model. More accurate 
measurements of the CMB, especially in the large multi- 
pole region, should give tighter constraints. Based on the 
data we cannot conclude that this model is ruled out. It 
would be necessary for formulate a specific model before 
more progress can be made concerning the question of a 
possible negative equation of state of dark matter. 

How strongly dependent are these results on our as- 
sumption of flatness and dark energy equation of state 
Wde = — 1? From Fig. [21 we can readily see that the main 
constraint on the r = model comes from the large scale 
structure data; since the CMB spectrum does not change 
much in the allowed parameter range, we may conclude 
that relaxing the flatness assumption would not make 
much difference on the constraints. The same is true for 
an equation of state different from — 1. From previous 
studies it is known that for Wde > — 1> structure growth 
is suppressed at small scales 0, 0, 0. But for the 
r = model, this can only make a small difference, since 
the growth suppression cannot ameliorate the strong de- 
viation from the LSS measurements at small scales as is 
readily apparent in Fig. 

The situation for the c^ d = case is different. Here, 
relaxing the flatness assumption and allowing for open or 
closed universes will lead to a weaker constraint on Wdm- 
The first peak position is sensitive to the geometry of the 
universe, but increasing or decreasing Wdm can in princi- 
ple shift this peak to be in agreement with the WMAP 
data, as may be seen in Fig. [31 We therefore expect 
also that the constraint on the total energy f2 tot will be 
weaker than in the standard ACDM case. Concerning 
the possibility that the equation of state of dark energy 
Wde > — 1 we may say that here, too, the constraints 
on the dark energy equation of state and Wdm will be 
less stringent. As mentioned above, the main impact of 
Wde > — 1 is through a suppression of structure growth at 
small scales, but this may be compensated by decreasing 
Wdm (see Fig. [3J|. It is therefore apparent that allowing 
for Wde =/= — 1 will result in weaker constraints on Wdm- 

Relaxing the flatness and Wde = — 1 will therefore have 
a negligible impact for the r = model but may lead 
to a significant relaxation of constraints for the (? ad = 
model. 




FIG. 4: One dimensional marginalized likelihoods of the MCMC simulation for the model parameters, using CMB, SNe la and 
FSS data. The results are for the T = (black, straight) and the c 2 ad = (red, dashed) model. We have plotted the results 
for a pure ACDM model with w = for comparison (blue, dotted). Note the large difference in constraints on w for the two 
models. Also, the F = model has nearly the same parameter distributions as a pure ACDM model. 
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